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Formulas  are  developed  for  accurate  evaluation  of  the  normalized 
Hallen-System  matrix  that  arises  while  solving  for  the  current  distribu- 
tion on  a wire  antenna  or  scatterer.  The  technique  developed  here  may 
be  generalized  to  other  similar  electromagnetic  problems.  The  thin 
wire  example  is  chosen  for  illustrative  purposes  and  also  for  comparison 
with  previously  known  representations  for  the  matrix  elements.  It  is 
observed  that  the  present  method  is  applicable  whenever  the  kernel  func- 
tions are  analytically  Fourier  (or  Laplace)  transformable. 
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I. 


Introduction 


The  Pocklington  form  of  the  integral  equation  for  the  current  dis- 
tribution on  a straight  thin  wire  (receiving  or  driven)  is  well  known  and 
is  given  by 


L -27r  -jk  R 

/ z')  - ad<^'dz'  =-jwe 

/ 47rR  o z 


(1.  1) 


With  reference  to  figure  1,  L and  a are  the  length  and  radius  of  the 
wire.  If  we  assume  that  the  current  density  is  azimuthally  symmetric, 
the  total  current  across  the  antenna  cross  section  at  location  z is  given 
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Using  equation  1. 2 in  equation  1.  1,  we  get 
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where  the  kernel  is  given  by 


K(z  - z') 
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Figure 


I.  Geometry  of  the  Problem  Shov;ing  the  Zoning 
of  the  Thin  Wire 
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If  we  treat  the  integral  in  equation  1.  3 as  a function  of  z and  solve  the 
differential  equation,  llallen's^  form  of  integral  equation  is  obtained  as 


I I(z')K(z  - z')dz' 
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= A sin  (k  z)  + B cos(k  z) 
o o 


sin[k^(z  - z')]dz' 


(1.5) 


Implicit  in  equation  1.  5 is  a harmonic  time  dependence  of  the  type 
exp(ju)t)  which  has  been  Fourier  transformed.  Quite  often,  e.  g. , tran- 
sient analysis,  it  is  useful  to  introduce  complex  frequency  s,  in  which 
case  equations  1. 2 and  1.  5 become 
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and 
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I(z')K(z  - z')dz'  = A sinh(7z)  + B cosh(7z) 
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E^*^^(z')  sinh[7(z  - z')]dz' 
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with 


where 


7 = s/c  = propagation  constant 

= characteristic  impedance  of  free  space 


(1.6) 
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Implicit  in  equations  1.  6 and  1.  7 is  a time  dependence  of  the  form  exp(st) 

and  the  tilda  signifies  a two-sided  Laplace  transform.  It  can  be  easily 

verified  that  equations  1.  3 and  1.  5 are  special  cases  of  equations  1. 6 and 

1.  7 by  setting  s = jtJ . Hccause  of  the  differential  operator  in  equation 

1.  6,  equation  1.  7 is  more  amenable  for  machine  solution,  which  is 

2 

usually  achieved  by  the  method  of  moments,  e.  g. , Harrington.  This 
method  converts  the  integral  equation  1,  7 into  a system  of  linear 
equations 


[Z  ] [I  ) = [f  (z)]  ; p,  q = 1.  2,  3,  • • • , (n  + 1)  1 

p.  q p p I 

or  ](I  ] = [f  (z)]  (1.10) 

p.  q ^ 

which  is  then  solved  for  the  unknown  current  distribution.  We  shall  call 

this  Z matrix  as  the  Ilallen-System  (or  simply  Il-S)  matrix.  This  matrix 

2 

is  to  be  distinguished  from  the  generalized  impedance  matrix  and  associ- 

g 

ated  eigen  impedances  , both  of  which  have  a more  physical  interpreta- 
tion. If  we  use  pulse  functions  for  expanding  the  unknown  current  and 
delta  functions  as  the  testing  functions,  the  elements  of  the  normalized 
H-S  matrix  are  given  by 
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where  R = (z  - z')  + (2a  sin  <l>'  /2)  and  p,  q = 1,  2 • • • , (n  + 1) . 

pip  J 

[I  ) in  equation  1.  10  is  a column  matrix  made  up  of  zone  currents 
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The  right  hand  side  of  equation  1.  10  is  given  by 
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sinh  [7(z  -z')]dz'  (1.12) 
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A and  B are  constants  to  be  determined  by  imposinjf  the  end  conditions 

I(z  = 0)  = I = 0 and  I(z  =L)  = I = 0. 

1 (n+1) 


For  purposes  of  this  paper,  we  are  concerned  with  the  matrix 
whose  elements  are  given  by  equation  1.  11,  in  which  Zp  and  z' 
are  respectively  the  observation  and  source  points,  and  A is  the  cell 
width.  It  is  observed  that  [Zt^i  ] is  a Toeplitz  symmetric  matrix. 

q 

However,  if  the  current  distribution  is  of  interest,  using  I =1,  > = 0 

1 (n+1) 

and  rearranging  equation  1.  10  becomes 
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where  [Zto  ] is  [Zto  ] with  its  first  and  last  columns  replaced  by 

^^p.  q ^^p,q  I j 

j =-sinh(7Zp)  and  Z^^  (n+1)  “ " • l^p  1 a column 

matrix  made  up  of  A,  I-,  I„  • • • • I , B,  and  the  elements  of  the  column 
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matrix  IV^]  are  given  by 


E^'^*^(z')  sinh[7(z 
z p 


z')]  dz' 


(1.  14) 


Now,  however. 


^^Np, 


is  not  a Toeplitz  symmetric  matrix. 


Returning  to  equation  1.11,  let  us  specialize  it  to  the  diagonal  ele- 
ment. With  a change  of  variable  x = (z'  - z ),  it  is  given  by 
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This  integral  is  seemingly  singular  (as  x and  (j)'  approach  0)  and  the  two 
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approximations  used  by  Harrington  and  Tesche  will  be  reproduced  below 
and  later  compared  with  a more  accurate  evaluation  used  in  this  paper. 
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Equation  1.  16  is  derived  from  1,  15  by  approximating 


R(<j)')  = x + (2a  sini^i'/2)  a:  x + a 


1/2 


1/2 


and  keeping  the  first  two  terms  in  the  expansion  of  exp(-  7R) , whereas 
equation  1.17  retains  the  dependence  and  also  keeps  only  the  first  two 
terms. 
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II.  Matrix  Elements  by  Transform  Method 

In  this  section  we  shall  develop  formulas  for  more  accurate  evalua- 
tion of  the  general  element  in  the  normalized  II- S matrix  and  later 
specialize  it  to  the  diagonal  element  and  compare  results  with  those  of 
equations  1.16  and  1.17. 


Consider  the  <f>'  integral  of  equation  1.11 
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This  has  been  shown  by  Ilallen  to  have  a Fourier  inverse  transform 
representation  given  by 


K(z  - z ')  = 


— ^ y*  j K^(aVpT7' jd§ 


(2.  1) 


Because  of  the  symmetry  property  (z^  and  z'  are  interchangeable) , 
only  the  cosine  part  of  the  integral  contributes  so  that 


~ 1 /*°° 

K(z  - z')  = — 5 / cos  l?(z  - z')]  I (u)  K (u)d§ 
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Using  this  result  in  equation  1.  11  and  performing  the  z'  integral,  a 
general  matrix  element  is  given  by 

~ \ S ^ - z )]  ^ - d^ 

N 2^00  "’pq 


p,  q TT  0 


(2.2) 
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If  z and  z are  the  mid  points  of  pth  and  qth  cells,  respectively  (see 

p q 

figure  1),  they  are  given  by 

z = (p  - 1)A  with  p = 1,2,3,  ••••  (n  + 1)) 

P (2.3) 

z = (q  - 1)A  with  q = 1,2,3,  ••••  (n  + D) 

q 

It  is  noted  that  the  first  and  last  cells  are  of  width  a/2  Jind  not  A which 
makes  the  end  points  z = 0 and  L to  be  the  centers  of  imaginary  cells  of 
width  A which  extend  beyond  the  physical  ends  of  the  antenna.  This  fact, 
if  and  when  necessary,  is  easily  accounted  by  dividing  the  first  and  last 
column  elements  by  a factor  of  2.  Substituting  equation  2.  3 into  equation 
2.  2 and  with  a change  of  variable  y = 5A/2 

p,  q TT  0 q.  p 

where  the  argument  v of  the  modified  Bessel  functions  is  given  by 


v 


2 

a 7 + 


For  adequate  representation  of  the  current,  if  we  choose  the  cell  width 
A » a and  « radian  wavelength,  then  equation  2.  4 is  quite  accurate.  The 
radian  wavelength  H may  be  computed  via  X = c/j  s|  . Equation  2.4  is 
now  ready  for  machine  integration  and  may  be  performed  for  real  y or 
one  can  treat  y as  a complex  variable  and  suitably  deform  the  contour 
on  the  real  axis  to  a contour  that  wraps  around  the  branch  cut.  The 
branch  points  are  at  y =lj7A/2.  We  shall  now  specialize  equation  2.4 
for  p = q and  obtain  the  normalized  diagonal  element  as 
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2 

TT 
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v = a 7 4 (4y' 
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III.  Numerical  Results 


To  perform  machine  integration  of  equation  2.  5,  it  is  useful  to 
investigate  the  singularities  and  the  analytical  structure  of  the  integrand. 
As  in  equation  2.  4,  the  integrand  in  2.  5 also  has  a pair  of  branch  points 
at  y =±j7A/2.  If  we  exclude  the  DC  frequency(T  = 0),  the  integration 
may  be  carried  out  on  the  positive  real  axis  of  the  complex  y plane.  The 
infinite  integral  of  equation  2,  5 may  be  broken  up  as  follows 


^DN 


I (v)K  (v)dy 
y o o 


(3.  1) 


If  we  choose  an  upper  limit  UL  which  permits  a large  argument  approxi- 
mation for  the  modified  Bessel  functions,  e.  g. , | v|  = 10  implies 
UL  = 5A/a  because  of  [yA  | « 1 . One  can  now  write 


I (v)  K (v)  dy  + 
y o o 


5A/a 


sin  y 

y 


(3.2) 


For  a suitable  change  of  variable,  it  is  possible  to  write  the  second  inte- 
gral over  a finite  range  of  integration.  Both  of  these  integrals  were 
carried  out  with  a 40-point  Gaussian  quadrature  routine.  Care  was  taken 
to  ensure  outgoing  wave  nature  by  choosing  the  proper  Bessel  function 
product,  depending  on  their  arguments.  The  ranges  of  integration  were 

continually  subdivided  until  a convergence  criterion  of  the  type,  the  mag- 

-3 

nitude  of  the  ratio  of  two  successive  answers  was  less  than  10  was  met. 

With  this  type  of  convergence,  it  was  seen  that  the  second  integral  was 

0 

insignificant  since  it  contributed  values  of  the  order  of  10  . Requiring 

a 3-figure  accuracy,  the  second  integral  was  neglected.  Numerical  com- 
putations were  made  for  cell  sizes  ranging  from  1 to  100  radii  for  the 
case  of  50  cells  per  wavelength.  In  order  to  be  able  to  compare  with 
previous  results,  the  computation  is  done  on  the  imaginary  axis  of  the 


-10- 


complex  s plane.  The  results  are  found  in  table  1 and  plotted  in  figure  2. 

2 

The  real  part  of  from  Harrington's  approximate  formula  is  a 

straight  line  in  this  figure  and  is  valid  only  for  A > 1 Oa,  The  real  part 

3 

from  Tesche's  formula  is  not  plotted  because  of  its  closeness  with  the 

values  obtained  by  the  transform  method.  This  can  be  seen  in  table  1. 

2 3 

Furthermore,  the  imaginary  parts  of  both  Harrington  and  Tesche  are 
eonstai^t  and  equal  to  each  other  because  of  similar  approximations  in- 
volved. It  is  seen  that  the  imaginary  part  from  the  transform  metliod,  as 
may  be  expected,  is  oscillatory  and  appears  to  have  a mean  value  equal  to 
the  constant  value  - 7A/(47r)  in  the  previous  two  methods.  It  must  how- 
ever be  pointed  out  that  both  of  the  previous  results  are  excellent  approxi- 
mations within  the  ranges  of  their  validity  and  the  transform  method, 
although  significantly  more  time  consuming  (10  times  or  more  in  several 
cases),  does  lead  to  improved  accuracy. 
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A/a 

Harrington 

Tesche 

Transform  Method 

Heal 

Imaginary 

Real 

Imaginary 

Real 

Imaginary 

1.0 

0.00000 

-. 01000 

. 09529 

-.01000 

. 09447 

-.00999 

2.0 

1 

.11032 

-.01000 

.15460 

-,  01000 

. 15472 

-. 00982 

3.  0 

. 17485 

-.01000 

. 20000 

-.  01000 

. 19974 

-. 01004 

4,  0 

. 22064 

-. 01000 

.23661 

-.01000 

. 23663 

-.00956 

5.  0 

. 25615 

-. 01000 

. 26709 

-.01000 

, 26660 

-, 00966 

6.  0 

. 28517 

-.01000 

. 29309 

-.01000 

. 29267 

-. 01028 

7.  0 

. 30970 

-, 01000 

, 31568 

-.01000 

. 31568 

-. 01051 

8.  0 

. 33095 

-, 01000 

. 33561 

-.01000 

. 33534 

-.01040 

9.0 

, 34970 

-. 01000 

. 35343 

-.01000 

. 35281 

-. 00999 

10.  0 

. 36647 

-. 01000 

, 36952 

-.  01000 

. 36946 

-, 00931 

20.  0 

. 47679 

-, 01000 

. 47757 

-.01000 

. 47691 

-.00956 

30.  0 

. 54132 

-. 01000 

, 54167 

-.01000 

. 54103 

-.01073 

40.  0 

. 58710 

-, 01000 

. 58730 

-.01000 

. 58630 

-.01015 

50.  0 

. 62262 

-. 01000 

. 62274 

-.01000 

. 62323 

-. 00834 

60.  0 

. 65164 

-.01000 

. 65172 

-.01000 

. 65053 

-. 01001 

70.  0 

. 67617 

-. 01000 

. 67623 

-.01000 

. 67598 

-. 01168 

80.  0 

. 69742 

-.01000 

, 69747 

-.01000 

: . 69848 

1 

-. 00782 

90.  0 

. 71617 

-. 01000 

. 71621 

-.01000 

. 71526 

-.00879 

100.0 

. 73294 

-. 01000 

, 73297 

-.01000 

. 73150 

-. 00977 

Table  1 

Diagonal  term  of  the  normalized  II alien- System  matrix  for  a 
thin  wire  as  a function  of  normalized  cell  width; 

X/A  = 50,  00  or  kA  = . 125664  . 
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Figure  2.  Plot  of  tbs  Diagonal  Term  of  the  Normalized  Hallen 
System  Matrix  as  a Function  of  Normalized  Cell  Width 
(A/a),  KA  = 0.125664  (or  50  Cells  Per  Wavelength) 


IV. 


Conclusions 


An  accurate  way  of  computing  the  impedance  matrix  elements  for  a 
thin  wire  structure  is  outlined.  This  procedure  is  based  on  being  able  to 
analytically  Fourier  or  Laplace  transform  the  kernel  function.  The  thin 
wire  structure  was  only  chosen  for  illustrative  purposes,  although  the 
method  is  applicable  if  and  whenever  the  kernel  is  analytically  transform- 
able. Another  example  where  this  method  is  usefully  employed  may  be 

found  in  the  treatment  of  three  dimensional  EM  scattering  from  a finitely 

4 5 

long  circular  cylinder  by  Kao.  ' 

Numerical  results  for  a representative  case  of  a thin  wire  structure 
are  compared  with  previously  available  results.. 
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